Coherent States Formulation of Polymer Field Theory 
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We introduce a stable and efficient complex Langevin (CL) scheme to enable the first numeri- 
cal simulations of the coherent-states (CS) formulation of polymer field theory. In contrast with 
Edwards' well known auxiliary-field (AF) framework, the CS formulation does not contain an em- 
bedded non-linear, non-local functional of the auxiliary fields, and the action of the field theory 
has a fully explicit, finite-order and semi-local polynomial character. In the context of a polymer 
solution model, we demonstrate that the new CS-CL dynamical scheme for sampling fluctuations 
in the space of coherent states yields results in good agreement with now-standard AF simulations. 
In addition, the new framework makes accessible the asymptotic ground-state limit of collections of 
infinitely long polymer chains with fluctuations included. The formalism is potentially applicable to 
a broad range of polymer architectures and may facilitate systematic generation of trial actions for 
use in coarse-graining and numerical renormalization-group studies. 



The standard auxiliary-field (AF) polymer field the- 
ory introduced by Edwards in the 1960s [l[, has proven 
to be a powerful tool for investigating the properties of 
mesoscopic models of inhomogeneous polymers at equi- 
librium and has been well developed on both analytical 
and numerical fronts over the past five decades. Apply- 
ing the mean-field approximation to Edwards' framework 
results in self-consistent field theory (SCFT), which has 
been the basis for exploring the self-assembly of a wide 
range of important heterogeneous polymer systems @B| ■ 
Moreover, field-theoretic simulations (FTS), which move 
beyond the confines of the mean-field approximation, of- 
fer a scalable and general computational methodology 
for studying fluctuation phenomena in Edwards' AF-type 
models, with numerical methods for circumventing the 
sign problem based on complex Langevin (CL) sampling 



The standard approach to formulating the Ed- 
wards AF framework involves the use of Hubbard- 
Stratonovich transformations to decouple segmental in- 
teractions through the introduction of one or more 
spatially- varying auxiliary potential fields 0"11|. The 
effective Hamiltonian resulting from this procedure is a 
nonlinear functional of the auxiliary fields, with the non- 
linearity embedded within a single-chain partition func- 
tion that enumerates the conformational states of non- 
interacting polymer chains subject to an external field. 
The non-linear and non-local character of such Hamilto- 
nians causes difficulties for advanced applications of poly- 
mer field theory, including systematic coarse graining and 
numerical renormalization group studies of a wide range 
of polymeric fluids, and investigations of rigid or semi- 
flexible liquid-crystalline polymers. 

In 1969, Edwards and Freed introduced an alternative 
formulation of polymer field theory based on the second- 
quantized description of the quantum many-body prob- 
lem, and adapted it to analytical studies of the statisti- 



cal thermodynamics of an assembly of crosslinked poly- 
mers 12j. However, the formalism has been largely ne- 



glected since that time, has not been explored in the 
context of conventional linear (non-crosslinked) polymer 
solutions and melts, and to our knowledge, has never 
been applied as a basis for numerical simulations. The 
distinguishing feature of this formalism is the absence 
of auxiliary external fields; rather, the only fields enter- 
ing the effective Hamiltonian are propagator-like coher- 
ent states (CS). Interestingly, the highly non-local char- 
acter of the auxiliary-field formalism is absent: the ac- 
tion in the CS formalism is a fully explicit, semi-local, 
fourth-order polynomial in the fields. One might hope 
that this aspect of the CS theory would facilitate simpler 
extension to branched and networked polymers, rigid or 
semi-flexible chains, and to coarse graining and renor- 
malization group studies. In this Letter, we revisit the 
Edwards-Freed CS formalism in the context of a stan- 
dard polymer solution model and introduce a stable and 
accurate stochastic dynamics scheme for conducting fully 
fluctuating simulations in the CS framework. 

We use the Edwards model of homopolymers in an 
implicit solvent (Model A in Ref. Q) as a platform to 
develop the algorithms necessary to conduct simulations 
in the CS formalism. The CS description can be de- 
rived from the auxiliary-field representation of the grand- 
canonical ensemble [3j . For a collection of homopolymers 
in implicit solvent, the grand partition function can be 
written 

SV .V.T>-± J 

n=0 J 

cxp^-i j dr J dr'wW^dr-r'Dw^j (1) 

where Q [w] is the single-chain partition function in the 
presence of the auxiliary field w(r), u(r) is the implicit 
pair potential acting between monomers in units of fcsT, 
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and z' is the chain activity. The partition function Q[w] 
can be expressed as a linear functional of a Green func- 
tion, Q[w] = / dr J dr' G (r, N; r', 0), where r and r' rep- 
resent the coordinates of the two polymer ends and N 
is the polymerization degree. By definition, G (r, s; r', s') 
is the inverse of the operator «5f = d s — ^-V 2 + iw ap- 
pearing in the standard diffusion equation for a poly- 
mer chain propagator Q ■ By recognizing that the Green 
function can be written as a correlator of fields p and <p, 
G (r, s; r', s') = — « (</? (r, s) <p (r', s')) , where causality is 
implied (non-zero for s > s'), and ip, tp appear in the ac- 
tion of a generating functional, Hq — — i J ds J dr ipJZ'tp, 
one may transform the sum over all powers of the Green 
function to 



^ w 



n=0 



JdrJ dr'G{r,N-r',0)Y 



JVpJVpe~ H " +VF I 



dr[-iip(r,0)+ip(r,N)] 



JVtpJVp e~ Ho 



(2) 



In the latter, the functional integrals, Dtp and T>tp, are 
understood to be over all possible domain-supported real- 
valued fields. Inserting Eq. ((2} into Eq. ([I} and complet- 
ing the Gaussian functional integral over w(r) leads to 



E(z',V,T) 



J Vtp J Vip exp (-H [<p, ip}) , 



(3) 



with effective Hamiltonian (action) 
H [tp, (f] = -i J ds J drp> (r, s) (d s - y V' 



p(r,s) 



+ - I dr I dr'p(r)u(\r-r'\)p(r') 
z' I dr [-ip (r, 0) + p (r, N)\ , 



where p is a monomer density operator given by 



p (r) = —i I ds 
Jo 



(p(r,s)(p(r, s). 



(4) 



(5) 



The terms in H [tp, tp] have an intuitive interpretation, 
with the first describing non-interacting polymers, the 
second term their pairwise excluded-volume interactions, 
and the third a source injection at the two chain ends 
necessary to generate the grand canonical ensemble of 
polymers. 

We note that the original Edwards formulation of this 
model employed delta function excluded-volume interac- 
tions for u(r). This choice, in which u is not finite on 
contact, leads to an ultraviolet (UV) divergent theory, in 
which averages of field-theoretic operators do not have 
a finite continuum limit 13 , [l4 1 . To regularize the the- 



ory, we instead adopt the repulsive Gaussian potential 
u(r) — [u /(87r 3 / 2 a 3 )] exp[— r 2 /(4a 2 )] with strength uq 
and range a. 



Developing viable numerical schemes for simulation of 
polymer field theories requires care in handling the sign 
problem that arises from complex actions. In this re- 
gard, the CS formalism is no different: straightforward 
integration over real CS fields, e.g. by Monte Carlo sam- 
pling methods, is prohibitively inefficient due to the sign 
problem. However, since the integrand is analytic in both 
fields, it is conceptually possible to deform the integra- 
tion paths to pass through complex- valued saddle points, 
which represent mean-field or SCFT solutions. Beyond 
the mean-field approximation, we have found that the 
complex Langevin (CL) dynamics schemefl, 0, [lj| is 
a useful technique for simulating complex-valued (AF) 
polymer field theories, since it is adaptive and does not 
require saddle points or constant phase contours to be 
computed in advance. While the above summarizes our 
experience with the AF formalism, we shall see that ad- 
ditional complications arise in devising a thermodynam- 
ically consistent CL scheme for the CS framework. 

For numerical convenience, we rescale s to [0,1], ex- 
press all spatial lengths in units of the free-polymer ra- 
dius of gyration R g = by/W/6, and absorb a factor of 

3 /2 

R g into both fields, resulting in the rescaled Hamilto- 



H [tp, <p] = -i J ds J dr tp (r, s) (d s - V 2 ) p (r, s) 
+ ^JdrJ dr'p(r)T(\r-r'\)p(r') 
-V*/"drH0(r,O) + ¥>(r,l)], (6) 



where B — uqN 2 / Rg is a dimensionless excluded volume 
parameter, z = z'R^ is a dimensionless chain activity, 
r(r) = [l/(87r 3 / 2 a 3 )] exp[— r 2 /(4a 2 )] is a reduced poten- 
tial function with range a = aj R g and unit volume inte- 
gral, and the dimensionless monomer density operator is 

Next we turn to discuss the analytic structure of the CS 
field theory and propose a dynamical relaxation scheme 
for numerical calculation of saddle points (SCFT) and 
a related CL dynamics to importance sample the func- 
tional integrals of the full field theory. It is evident 
from Eqs. (J3HSI) that any saddle point solutions, with 
associated real-valued observables such as the density 
operator, must lie off the integration path; for exam- 
ple with one CS field becoming purely imaginary. We 
therefore promote the CS fields to be complex. The ac- 
tion is everywhere analytic, so deformation of the inte- 
gration path poses no immediate problems. A simple 
complex dynamics that has saddle points as its steady 
state is the "diagonal descent" scheme: d t p = —SH/8p>, 
dtp = —5H/Sip. These dynamics indeed recover phys- 
ically realistic saddle-point solutions; specifically, solu- 
tions where tp is purely imaginary and tp is real. In the 
case of periodic boundary conditions, it is straightfor- 
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ward to show that the only saddle point/SCFT solution 
of the present polymer solution model is spatially ho- 
mogeneous and given by (p(s) = </z exp(— Bpms)Q(s), 
<p(s) — iy/zex.p[—Bp M (l — s)]0(l — s), where 0(s) 
is the Heaviside function and the mean-field reduced 
monomer density pu is given as the solution of pu — 
z exp(— Bpm)- This mean-field solution is of course iden- 
tical to SCFT analysis of the familiar AF representation 
of the model. In deriving these results, or implementing 
numerical field relaxation (or CL below) in less trivial 
situations, we note that the source terms in Eq. ^ pro- 
duce delta function forcing at the two chain ends; these 
are best treated as (time-independent) boundary condi- 
tions, namely <p(r,0) — yfz and <p(r, 1) = i\fz. 

Unfortunately the diagonal descent scheme is numeri- 
cally ill-behaved, exhibiting marginal stability and oscil- 
latory trajectories to saddle points. It is therefore un- 
suitable for obtaining SCFT solutions and as a basis for 
stochastic CL simulations. However, a stable scheme is 
obtained by applying a non-Hcrmitian mobility matrix 
to the field relaxation equations: 



dt 



i 

1 



( SH/Sp \ 
V SH/S0 J 



(7) 



where we have also added complex- valued, Gaussian, and 
white (decorrelated both spatially and temporally) noise 
sources p, and p to be used in the case of CL dynamics. 
This matrix has the effect of decoupling the dynamics for 
ip and (p to linear order, while also implicitly including 
a Wick rotation for stable saddle-point relaxation. It is 
readily verified that without the noise terms, Eq. ([7]) is 
a robust scheme for relaxing SCFT solutions of the CS 
field theory. 

Beyond mean field theory, the noise terms fi and (i 
must be carefully constructed in order to produce a CL 
dynamics scheme that is thermodynamically consistent; 
namely that has a steady state probability distribution 
of field states consistent with the complex Boltzmann 
weight of the partition function in Eq. ([3]). As complex 
Langevin theoryjl, [3] has not been developed for CL 
equations with non-Hermitian mobilities, we provide a 
proof in the Online Supplemental Information that ther- 
modynamic consistency imposes six restrictions on the 
second moments of the real (R) and imaginary (I) com- 
ponents of the noise terms (the first moments vanish) : 
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< fi% > - < A*? > = < Afl > - < A? > 

< m.r£.r > = < M/Aj > 

< (LtRjUj > = < fl R fll > = 







(8) 



and 



< /i fl (r,s,i)/2i(r', s',t') > + < /i/(r, s, t)p, R (r', s', t') > 
= 26(r-r')S(s-s')6(t-t').(9) 

These conditions are insufficient to uniquely specify the 
noise statistics for a CL model. While all choices con- 



sistent with the above constraints should yield identi- 
cal time average properties, we expect that trajectory- 
specific properties such as population variance and cor- 
relation times can vary significantly. Among a vari- 
ety of schemes investigated, our best performing CL 
model uses the method described in [17| to construct 
thermodynamically consistent, correlated complex noises 
from four real and independent Gaussian white noises, 
iji, i — 1,2,3,4, with vanishing first moments and 
second moments given by (rji (r, s;t) rjj (r', s', t')) — 
26 (r — r') 5 (s — s') S (t — t') 6ij. The complex noise 
components are p, R = r)x, \ij = i] 2 , fiR = (r/ 2 + n/3^) /2, 
and /tj = (jii + VSn^ /2. Combining these noise com- 
ponents with Eq. ([7]). we arrive at the following CL dy- 
namics for conducting field-theoretic simulations in the 
CS framework: 



.SH 

dip 

n a .SH rj 2 + V3n 3 . 
dt(fl = + 2 +l 



(10) 



We now demonstrate that a numerical implementation 
of this CL dynamics scheme has at steady state time 
averages that correspond to the correct ensemble aver- 
ages of the model defined by Eqs. For this pur- 
pose, plane wave spectral collocation was used in a cubic 
simulation cell subject to periodic boundary conditions. 
Second-order finite differences were applied in the con- 
tour coordinate s and a weak first-order semi-implicit 
integration method [15( was used to time step the dis- 
cretized, stochastic CL equations. All calculations were 
conducted in a cell of side length L = 3.2 R g with B = 1, 
Ax = a = 0.1 R g , As = 0.01, and At = 0.01. Since exact 
results are not available for the fully fluctuating theory, 
we also conducted simulations using the conventional AF 
representation of the model, with identical discretization 
procedures and parameters (14j. 

Fig. Q] shows the fluctuation contribution to the spa- 
tially averaged polymer density, (pr), relative to the 
mean-field value pm, for various values of the chain ac- 
tivity, z. The cell size of L = 3.2 R g was chosen to elim- 
inate finite-size bias from truncating below the solution 
correlation length, as shown in the inset of Fig. [T] The 



cutoff, 



0.1 R„, was similarly chosen as a microscopic 



scale to regularize the UV divergence in the theory with- 
out changing properties on relevant mesoscopic length 
scales, i.e. R g or the correlation length, while Ax and 
As were selected to provide a well-resolved spatial and 
contour sampling of such a regularized field theory. Fi- 
nally, At = 0.01 was chosen on the basis of a convergence 
study (see Online Supplemental Information) to be small 
enough to achieve stable and accurate field trajectories. 
The general agreement and statistical consistency seen 
between simulation results for the AF and CS represen- 
tations of the model suggest that the CS framework is 
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a viable platform for conducting field-theoretic polymer 
simulations. 
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FIG. 1. (Color online). Comparison of the values of the real 
part of the spatially and temporally averaged polymer den- 
sity, < pr >, which have been subtracted from their corre- 
sponding saddle-point values, pM, for different values of the 
activity z. Statistical consistency is seen between the results 
from simulations based on Edwards' auxiliary field represen- 
tation and the new coherent states formalism. Calculations 
were conducted for 5 = 1, Aa; = a = 0.1 R g , L = 3.2 R g , 
As = 0.01, and At = 0.01. The error bars displayed are the 
stochastic error of the mean. The inset shows suppression of 
finite-size bias on pM— < pn > by varying L from 1.4 R g to 
5.0 R g with all other parameters fixed. 



Beyond the fact that the CS description has a simpler 
and more local action than the AF model representa- 
tion, which could be advantageous in numerical coarse- 
graining or RG studies, the CS formalism also facili- 
tates the study of fluctuation phenomena in collections 
of asymptotically long polymer chains. A good example 
is a semidilute polymer solution, where the structure and 
thermodynamics are dictated by a molecular weight in- 
dependent correlation length £, while the overall polymer 
size, i.e R g , can be orders of magnitude larger. The CS 
framework provides a natural vehicle for simulating such 
situations without the explicit need to resolve both the £ 
and R g scales. 

To investigate collections of long polymers, we impose 
a ground state dominance (GSD) approximation within 
the CS framework. This proceeds by expanding the CS 
fields in the complete basis of eigenstates of Jz? 



ip(r,s) = ^=y^e~ 

i=5>^(r) 



ip(r,s) 



(11) 



K 



where fix is Kth eigenvalue and ipK (r) is the correspond- 
ing eigenstate. In the GSD approximation, which is valid 
in the limit of very large N, we retain only the lowest 
eigenstate, \i m : 



1 



Inserting Eq. (172]) into Eq. (J4j) yields 



H[^,fl>] = -i I dr^j (r) ( - f i m - -V 2 ) V (r) 



+ ljdrj dr'p(r)u(\r-r'\))p(r') 

-#(r) + e^^r)! , (13) 



s ip(r), <p(r,s) 



N 



V>(r)(12) 



where the monomer density is now p(r) = —iip(r)ip(r). 
Taking N — > oo, it is clear that the source term in 
Eq. p^|) disappears, resulting in a GSD action, 

H$, — i j dv 'tp (r) (-fi m - 1 V 2 ) V (r) 

+ T / dr J dr ' p (r) r (|r " r ' l} p (r ' } ' (14) 

where all lengths have been rescaled by the statistical 
segment length b, and B' = uo/b 3 is a dimensionless 
cxcludcd-volume interaction parameter. Note that \i m 
takes the role of a monomer chemical potential, which is 
the only relevant quantity in a system of asymptotically 
long chains. 

Using the CL framework described above, it is straight- 
forward to conduct simulations of the GSD approximated 
CS theory - simulations that would be difficult or impos- 
sible in the AF framework. In Fig. [21 we report results of 
mean-field and fluctuating (CL) simulations within the 
GSD approximation for the monomer density at various 
monomer chemical potentials, conducted in a cubic sim- 
ulation cell of side L = 326 with B' = 1, Ax = 1, a = 5b, 
and At = 0.01. The mean- field results are different from 
the CL results at small values of fi m , indicating that for 
low-density systems, fluctuation effects can not be ig- 
nored even in the ground-state, long-chain limit. How- 
ever, as discussed in the Online Supplemental Informa- 
tion, the GSD CS model possesses UV divergences even 
with finite-contact excluded-volume interactions, which 
is manifest as a sensitivity to Ax in CL simulations of 
the model. 

To conclude, we have introduced the first numerical 
study of a coherent states polymer field theory. For this 
purpose, a complex Langevin scheme was devised that 
circumvents the sign problem, and is both stable and ef- 
ficient. Future developments will explore the use of the 
CS framework for conducting systematic coarse-graining 
and numerical renormalization group studies [rsl]. as the 
locality and finite polynomial order of the theory should 
simplify the identification of suitable basis functionals 
to be used in trial coarse-grained actions. The coherent 
states representation may also provide advantages in con- 
ducting field-based simulations of assemblies of branched 
and networked polymers or interacting semiflexible poly- 
mer chains. 

We acknowledge support of the National Science Foun- 
dation under grant DMR-CMMT-1160895 and computer 
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FIG. 2. (Color online). Comparison of the average values of 
the monomer density, {pit), between fluctuating and mean- 
field simulations in the ground state limit for various values 
of the monomer chemical potential, /i m . The error of the 
mean for the CL data is smaller than the symbol size. For all 
calculations B' = 1, Ax = 1 b, a = 5 6, and At = 0.01 with 
simulation cell size L = 32 b. 
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